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ABSTRACT 

We discuss the recent progress in our understanding of neutrino flavour con- 
version in core-collapse supernovae and focus on the effects coming from the 
neutrino-neutrino interaction. The latter has been shown to engender new phe- 
nomena, modifying significantly the neutrino fluxes and supernova observations. 
In particular, we make the link between the spectral split and the magnetic 
resonance phenomenon. Moreover, recent investigations have demonstrated the 
existence of leptonic CP violation effects on the supernova neutrino fluxes. We 
summarize these findings. 

1. Major progress ongoing... 

In the last decade simulations of core-collapse supernova explosions have reached 
an impressive complexity with the extension to two- and to three- dimensions, the 
inclusion of convection and of detailed neutrino transport. This has revealed the 
role played by the Standing- Accretion-Shock Instability (SASI mode). At present 
simulations for different progenitor masses show indications that an explosion occurs 
due to a com bination of convection, shock-wave instabilities and neutrino energy 
depositiorP^ Assessing the role of magnetic fields and of rotation needs further 
investigation as well. The longstanding question of how the life of massive stars ends 
is likely to be answered in the coming decade. 

Our interest in these explosive environments also resides in the gigantic amount 
of neutrinos taking away most of the gravitational energy in a few seconds. On one 
hand measuring the neutrino luminosity curve will give us an image of the different 
phases of the explosions, from the neutronization burst to the cooling of the proto- 
neutron star, or to black-hole formation. On the other hand, we hope to get clues to 
key neutrino properties that remain unknown, such as the neutrino mass hierarchy. 

Besides the interest for supernova simulations and for unravelling future observa- 
tions, the study of neutrino flavour conversion in supernovae is of intrinsic theoretical 
value. In fact, when we think of how neutrinos change their flavour while travelling in 
matter, the process that has become a reference in our minds is the Mikheev-Smirnov- 
Wolfenstein (MSW) effect. It is now established by the beautiful results of the SNO©, 
KAMLAND® and recently Borexino^ experiments. The MSW effect consists in an 
enhanced neutrino flavour conversion that originates from the fact that neutrinos un- 
dergo an adiabatic resonant conversion in the sun, due to their coherent scattering 
upon electrons. It is at the origin of the solar neutrino deficit. However, in the last 



years, it has become apparent that flavour conversion in core-collapse supernovae in- 
volves unexpectedly more complex phenomena. This complexity has emerged after 
realizing calculations of the neutrino propagation including the neutrino-neutrino in- 
teraction contribution to the Hamiltonian, using temporally evolving matter density 
profiles with shocks and with matter density fluctuations. As a consequence, predic- 
tions of core-collapse supernova neutrino fluxes have become demanding. A realistic 
calculation nowadays tipically asks for the solution of a large number of stiff, coupled, 
non-linear equations. On the other hand, simplified models are often employed to pin 
down the underlying mechanisms observed in numerical simulations. While major 
progress has been achieved, many questions remain open and surprises might still 
well be around the corner. 

The regions where neutrino mixings plays a role is far out the neutrinosphere, 
where one can safely neglect incoherent scattering. The inclusion of such contribu- 
tions, that go like G\ (Gp being the Fermi constant), would require the solution of 
the Boltzmann equation, as done in implementing the neutrino transport in super- 
nova models. Such simulations provide us with the "primary" neutrino spectra at 
the neutrinosphere (the region where neutrinos start free-streaming). At present the 
primary neutrino fluxes turn out to be well described by (pinched) Fermi-Dirac or 
power-law distributions, characterized by the neutrino average energies and the pinch- 
ing. For these two parameters, still a range of values is allowed, which constitutes 
a supplementary source of uncertainties in predicting the supernova neutrino fluxes 
that might be observed on Earth. In the future the achievement of succesfull super- 
nova explosions will certainly shrink the current range. However, it is also important 
to search for (combination of) observables, in supernova observatories, that can help 
pinning down the primary neutrino fluxes, in spite of the ensemble of possible neutrino 
fluxes, associated with the unknown neutrino properties and with the different flavour 
conversion phenomena that can occur outside the neutrinosphere. Indeed, extracting 
information on such primary neutrino fluxes, from a future supernova signal in an 
observatory, would offer an important test of supernova simulations. To conclude, 
while at present it does not appear as compelling to implement the neutrino mixings 
in the Boltzmann treatment, merging the two descriptions could well be one of our 
goals in the very long run. 



2. A lot of work already done... 

At the neutrinosphere, the neutrino density is so large that the neutrino interaction 
with neutrinos becomes significant. It was first pointed out by Pantaleone® that the 
presence of the vv interaction introduces an off-diagonal refraction index. The first 
numerical calculations showing that new phen ome na might arise were performed by 
Samuel, while only recently the results of RefP^ have triggered intense theoretical 
activity. The flavour content of a neutrino emitted at the neutrinosphere might have 




Figure 1: The neutrino-neutrino interaction effects upon the neutrino evolution in a core-collapse 
supernova : Diagonal matrix elements of the neutrino Hamiltonian within two-flavors neutrino 
evolution and in the 'matter' basis (fci + Qi)/6kvi (upper), (&2 + Q^l^n (middle) and their 
difference 1 + SQi2/Ski2 (lower figures). The curves correspond to a 5 (black) and 10 MeV (grey) 
neutrino energy for neutrinos (left) and anti-neutrinos (right figures). The black (grey) lines show 
the average for a 5 (10) MeV neutrino. The calculations include the vacuum mixing, the coupling 
to matter and the neutrino-neutrino interaction. The abrupth change of the matrix elements is due 
to a rapid growth of the matter phase associated with the onset of bipolar oscillationi^ . 



changed at the point where it interacts with another neutrino. Implementing the vv 
interaction introduces non-linearity into the problem. As done for the MSW effect, 
in the current literature, the neutrino-neutrino effective Hamiltonian is based upon 
the mean-field approximation. Note however that a recent workP^ :) has shown the 
exact solvability of the many-body Hamiltonian for the case of two-neutrino flavours 
(in the absence of matter and within the single-angle approximation). 

Currently three flavour conversion regimes have been identified. These are called 
the synchronization region, the bipolar oscillations and the spectral split (see the 
revie and references therein). Nearby the neutrinosphere the neutrino- neutrino 
interaction term is so large that each neutrino polarization vector essentially precesses 
around the corresponding effective magnetic field in a similar way and flavour con- 
version is frozen. The onset of bipolar oscillations has been studied in several works 
see e.e; .[121). In RefPJ the neutrino evolution in presence of the neutrino-neutrino 
interaction term has been investigated in the 'matter' basis that instantaneously di- 
agonalizes the Hamiltonian with matter (and the vv interaction) term(s). This has 
shown that the onset of the bipolar oscillations is triggered by a rapid growth of 
the matter phase (and its derivative). (Note that in the MSW case only the evo- 
lution of the matter angle matters.) The matter phase indeed arises because the 
neutrino-neutrino interaction gives complex contributions to the off-diagonal terms 
of the Hamiltonian in the flavour basis. We have found that, it is the growth of the 




Figure 2: Polarization vector formalism: The figure depicts the effective magnetic field, its z- 
component and the component lying in the (XOY) plane with the corresponding frequencies. The 
neutrino evolution can be seen as the precession of the neutrino polarization vector P around an 
effective magnetic field. Within this formalism one can link the spectral split to the magnetic 
resonance phenomenorP^ . 



matter phase that triggers flavour conversion since it renders the off-diagonal con- 
tribution of the Hamiltonian, in the matter basis, comparable to the difference of 
the diagonal terms. In the last phase of the evolution, depending on their energies, 
electron neutrinos and anti-neutrinos can completely convert into muon and tau neu- 
trinos implying a complete swapping of the fluxes. This is called the spectral split 
phenomenon and the energy at which such a swapping occurs is called the spectral 
split energy. Note that the vv interaction effects are independent of the value of the 
third neutrino mixing angle, provided that #13 is strictly non-zero, while their final 
impact on the neutrino fluxes depend upon neutrino properties and supernova aspects 
(e.g. the neutrino mass hierarchy, the initial flux ratios at the neutrinosphere and 
how dense matter is). 

Using simplified models, several works have brought very valuable insights into 
the collective effects engendered by the neutrino self-interaction and have given a 
vision of what actually is seen in the full numerical calculations . The approximations 
often made are: (i) neglecting the matter term; (ii) using a constant coefficient for the 
neutrino self-interaction term or with a simplified time dependence; or (iii) employing 
a single neutrino or anti-neutrino energy or neutrino box spectra. It was first argued 
m RefP, that the bipolar oscillations could be understood as a flavour pendulum. 
RefP 

has gone further and shown that if a neutrino-antineutrino asymmetry is 
present, as is case for the supernova neutrino fluxes, the evolution can be described 
in terms of a gyroscopic pendulum. The polarization vectors undergo both precession 
and nutation. A detailed study based on the gyroscopic pendulum is furnished by 



RefP^ 1 analyzing both the single-angle and the multi-angle cases. 

Most of the available calculations of neutrino propagation in supernovae including 
the neutrino-neutrino interaction assume the so-called "single-angle" approximation 



where the flavour evolution history is trajectory independent; while "multi-angle" 
calculations consider the flavour history along different trajectories and different in- 
teraction angles are considered. While the former has been shown to catch well 
qualitatively and quantitavively many features of the multi-angle calculations 
the two calculations reveal differences that can be important. In particular in multi- 
angle calculations collective effects present flavour decoherence, as discussed e.g. in 
while, when the matter density exceeds (is comparable) the neutrino density, col- 
lective effects can be strongly suppressed (or be affected by multi-angle decoherence). 
A recent investigation^^ 1 using matter density profiles from realistic ID supernova 
simulations has shown that dense matter can suppress collective effects during the 
accretion phase. It is clear that the interplay between the neutrino-neutrino interac- 
tion and dense matter needs further investigation. In RefP 

the authors have also 

shown that the location where bipolar oscillations start varies if a single-angle or a 
multi-angle calculation is performed. The implications of different onset locations of 
the bipolar oscillations on the r-process has been further investigated in where 
it was shown that multi-angle versus single-angle calculations can produce different 
r-process abundances. Note that the importance of the neutrino-neutrino interaction 
for the nucleosynthesis of heavy elements was identified in an early work^. Ref.^Q 
represents the first calculation where the vv interaction and shock wave effects have 
been put together in a consistent way. An interesting signature has been found for 
the positron time signal in a detector using inverse beta decay. When the shock wave 
passes through the MSW region a dip (bump) appears midway through the time 
signal depending on the positron energy, if the hierarchy is inverted and #13 is large. 

3. The neutrino spectral split as a magnetic resonance phenomenon 

At present, for the spectral split several features have been clarified. 

RefP 

has considered the spectral split a pure precession solution, which has been further 
extended in RefP3 

where an adiabatic solution in the comoving frame is postulated. 
This allows to predict in particular the split energy using lepton-number conservation 
which is an exact invariant of the two-neutrino flavour Hamiltonian (if matter is 
neglected), as first shown in^^Sl. In three flavours the number of invariants increases 
while the identification of the invariants for the many-body Hamiltonian is done in 
^21. RefJ23J has shown that multiple splits can appear depending on the neutrino flux 
ratios at the neutrinosphere. A global understanding of the occurrence of the splits 
is not yet achieved. 

Recently in Ref.ES we have been pointing out that the spectral split phenomenon 
is indeed a magnetic resonance phenomenon. Let us consider a spin precessing around 
a constant magnetic field Bqz, which can be taken along the z-axis, with a precession 
frequency given by u . The spin is also subject to a second magnetic field located in 
the (xOy) plane, which sinusoidally varies in time with u frequency. Such a magnetic 
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Figure 3: Three dimensional contour plot of the electron neutrino survival probability (in the flavour 
basis) as a function the neutrino energy and distance within the supernova. The middle of the 
white band region (color yellow-green region) corresponds to a survival probability of 0.5 so that 
all neutrino energies below the split energy of 7.6 MeV do not undergo any spectral swap, while all 
energies above the split energy undergo a spectral swap. The blac k cu rve superposed shows when 
the magnetic resonance condition in the matter basis is met : Co = (J2SJ . 



field acts as a perturbation that can however significantly impact the dynamics of 
the system depending on the relation between Au> = u — u>q and ui, the precession 
frequency around the planar field. The effective magnetic fields in the XY-p\ane, 
B± are associated with the off-diagonal components; while Bz corresponds to the 
difference of the diagonal matrix elements of the neutrino Hamiltonian (Figure 2). 
In RefJ23 

we have performed a full numerical calculation of the precession fre- 
quencies in the case of two-neutrino flavours whose evolution is determined by the 
Hamiltonian including neutrino mixings, their coupling to a realistic supernova den- 
sity profile and the neutrino-neutrino interaction calculated self-consistently (within 
the single-angle approximation). The magnetic resonance conditions that need to be 
fulfilled, for an inversion of the neutrino polarization vectors, ar Aw = oj — uq = 
and Au/ui <C 1. (within the flavour basis); or oj = and u;/u)i <C 1 in the comoving 
frame. Here Co (cDi) is the precession frequency of the polarization vector around the 
Z-(planar-) component of the magnetic field in the matter basis (Figure 2). We have 
identified the comoving frame given by the numerical average of the (fast varying) 
magnetic field associated with the Hamiltonian in the matter basis. We have founco 
that the (anti-)neutrino whose energies fulfill the magnetic resonance condition are 
the same that undergo the spectral split, at the same location in the supernova for 

a Note that our interpretation of the sp ectral split in terms of a mag netic resonance phenomenon 

lOOIO^l 11 k i 

is in agreement with previous work a * ' and the subsequent work^-W . 




Figure 4: Connection of the spectral split phenomenon with the magnetic resonance phenomenon 
for neutrinos: The region within whited bands (colored blue) shows the neutrino energies and the 
supernova location for which the resonance conditions Co rj (black curve) and Co/Coi <C 1 (region 
within white bands) - close to resonance - are both fulfilled (see text). When the resonance criteria 
are fulfilled, the flipping of the neutrino flavour polarization vector occurJ^^ . 



which the spectral split occurs (Figures E] and 4). 

4. Leptonic CP violation effects in core-collapse supernovae 

The search for a non-zero Dirac phase will be addressed either with upgrades of the 
T2K and NOz/A e xperim ents or with next generation acceleration facilities: super- 
beams, beta-beam anc [ neutrino factories. While the discovery of leptonic CP 
violation can help explaining the neutrino-antineutrino asymmetry in the Universe, 
it is an interesting question to ask if a non-zero Dirac phase can have an impact 
in other environments. An example is furnished by the neutrino fluxes in a core- 
collapse supernova. Only recently it has been demonstrated^! that such effects 
indeed exist and can arise, e.g. due to either radiative corrections or to physics 
beyond the standard model like flavour changing neutral currents. More generally, 
any physics that can break a factorization condition of the HamiltoniarPSJ can induce 
CP violation effects on the supernova neutrino fluxes. This result, first obtained in 
the framework of the MSW effect, has been shown to be valid also in presence of 
the neutrino- neutrino interaction^ 1 . Numerical calculations show that a non-zero 
Dirac phase induces modifications of the supernova neutrino fluxed^ 1 at the levef] of 
5-10%. On the other hand, the impact on the electron fractional - a key parameter 
for r-process nucleosynthesis - is loweiEl than 0.1 %. While it is now clear that effects 
coming from a non-zero Dirac phase exist in core-collapse supernovae, determining if 
they can have a phenomenological impact requires futher investigations. At present, 



b This result is based on a calculation which includes the vv interaction. 

c This result is based on a calculation which does not include the vv interaction. 




Figure 5: One- and two- neutron emission rates in the HALO phase-II observatory for a galactic 
explosion. The different curves show predictions taking into account the current uncertainty on the 
neutrino fluxes at the neutrinosphere from supernova simulations and the unknown v properties. 
The calculations include the vv interaction, coupling to matter and decoherence (see Ref. f or 
details). This figure is an example of how having detectors with different energy thresholds, or 
combining detectors based on different technologies, will allow to exclude possible solutions. 



the effect of a non-zero CP violating phase on the primordial abundance of light 
elements is also being investigated^. 

5. Remarks on observational aspects 

Unravelling the information encoded in supernova observations represent a chal- 
lenging task and what we will be able to extract will rely upon the knowledge by 
the time observations are made. On the neutrino physics side, this will depend upon 
the meausurement of the third neutrino mixing angle and possibly the hierarchy 
while other non-standard aspects can have an impact, such as the existence of sterile 
neutrinos or of non-standard interactions. Although degeneracies in the possible in- 
terpretations of observations will certainly be present, it will be useful to exploit the 
complementarity among the time and energy information associated with the differ- 
ent neutrino flavours emitted in a supernova explosion, the detection of the diffuse 
supernova neutrino background and the nucleosynthesis processes where neutrinos 
play a role {y process, r— process, i/p-process). 

The availability of a network of core-collapse supernova observatories based on 
several technologies - having a sensitivity both to electron neutrinos, electron anti- 
neutrinos as well as neutrinos of all flavours is clearly necessary -(scintillators, liquid 

q -i 

argon, water Cherenkov, see e.g. Ref. 1 ^). Different energy thresholds are clearly 
an advantage since the detectors are then sensitive to different parts of the neutrino 
energy spectra that might reflect the neutrino flavour conversion phenomena occur- 
ring in the supernova. They can present e.g. multiple energy splits due to the vv 



interaction that e.g. a lead-based detector like HALO, currently under construction 
at SNOLAB, can be sensitive to. Figure 5 shows the one- and two-neutron event rates 
expected in HALO phase-II for a galactic supernova explosion. The curves take into 
account the current uncertainties on the primary neutrino fluxes, the vv interaction 
and MSW effects as well as decoherence. This gives an example of how, depending 
on the neutrino average energies, one should be able to discriminate among either 
different primary fluxes or ^-properties (see Ref. 1 ^ for details). 

In the perspective of optimizing the information from future observations a precise 
calibration of neutrino-nucleus interaction cross sections is necessary for nuclei that 
are employed as targets. Low energy neutrino scattering experiments would have 
several applications, from nuclear structure studies to precise tests of the Standard 



physics program at intense spallation sources or realized at a low energy beta-beam 



6. What is ahead? 

The domain of neutrinos and core-collapse supernovae is undergoing major progress. 
While important achievements have been made, serious work is still needed. The in- 
crease in complexity has triggered fascinating theoretical questions and demanding 
numerical calculations that will be necessary to interpret future observations. Many 
important open questions remain that will steadily increase the complexity in the- 
oretical predictions. Indeed further efforts are required to achieve a comprehensive 
understanding of neutrino flavour conversion phenomena in presence of the neutrino- 
neutrino interaction and their interplay with other dynamical features. In very long 
run, future calculations will be based upon the multi-angle approach and include 
shock-wave, turbulence, get rid of the spherical symmetry assumption and use real- 
istic matter density profiles from supernova simulations... Tackling these theoretical 
issues as well as setting up an observational strategy for the network of supernova 
observatories represents a challenge for the future. 

7. References 



1) A. Marek and H. T. Janka, Astrophys. J. 694 (2009) 664 ! arXiv:0708.3372l 

2) S. W. Bruenn et al. larXiv: 1002 .49141 

3) Q. R. Ahmad et al, Phys. Rev. Lett. 89 (2002) 0113 01 |nucl-ex/0204008[ 

4) K. Eguchi et al, Phys. Rev. Lett. 90 (2003) 021802 |hj>ex/0212021[ ~ 

5) G. Bellini et al, Phys. Rev. D82 (2010) 033006 larXTv:0808.2868l 

6) J. T. Pantaleone, Phys. Lett. B 287 (1992) 128. 

7) S. Samuel, Phys. Rev. D 48 (1993) 1462. 



Model 




Such experiments could be part of the 




H. Duan, G. M. Fuller and Y. Z. Qian, Phys. Rev. D74 (2006) 123004 , 
|astro-ph/0511275| 

9) H. Duan, G. M. Fuller, J. Carlson and Y. Z. Qian, Phys. Rev. D74 (2006) 
05014, astro-ph/06066T6l 

10) Y. Pehlivan, A. B. Balantekin, T. Kajino, T. Yoshida, larXiv: 1105.11821 

11) H. Duan, G. M. Fuller, Y. -Z. Qian, Ann. Rev. Nucl. Part. Sci. 60 (2010) 569, 
larXiv:1001.2799l 

12) S. Hannestad, G. G. Raffelt, G. Sigl and Y. Y. Y. Wong, Phys. Rev. D7 4 
(2006) 105010 [Erratum-ibid. D76 (2007) 029901], ja?Xiv:astro-ph/0608695| 

13) S. Galais, J. Kneller, C. Volpe, larXiv:! 102. llTTl 

14) M. -R. Wu, Y. -Z. Qian, Phys. Rev D84 (2011) 045009. la7Xiv:1105.2068l 

15) G. L. Fogli, E. Lisi, A. Marrone and A. Mirizzi, JCAP 0712, 010 (2007), 
larXiv:0707. 19981 

16) A. Esteban-Pretel, S. Pastor, R. Tomas, G. G. Raffelt and G. Sigl, Phys. Rev 
D76 (2007) 125018. la7XTv:0706.24981 

17) S. Chakraborty et al.„ Phys. Rev. D84, 025002 (2011). la?XTv: 1105. 11301 

18) H. Duan, A. Friedland, Phys. Rev. Lett. 106, 091101 (2011). larXi"v:1006.2359l 

19) H. Duan, A. Friedland, G. C. McLaughlin, R. Surman, J. Phys. G38 (2011) 
035201 . larXiv: 1012.05321 

20) A. B. Balantekin and H. Yuksel, New J. Phys. 7 (2005) 51, |astro-ph/0411159| 

21) J. Gava, J. Kneller, C. Volpe and G. C. McLaughlin, Phys. Rev. Lett. 103 
(2009) 071101. larXTv:0902.03171 

22) H. Duan, G. M. Fuller, J. Carlson, Y. -Z. Qian, Phys. Rev. D 75 (2007) 
125005. |astro-ph/070"3776|. 



23) G. G. Raffelt and A. Y. Smirnov, Phys. Rev. D76 125008 (2007), 
larXiv:0709.46411 

24) B. Dasgupta, A. Dighe, G. G. Raffelt and A. Y. Smirnov, Phys. Rev. Lett. 
103 (2009) 051105. lirXlv: 0904 .35421 

25) S. Galais, C. Volpe, larXiv: 1103. 53 02, appear in Physical Review D. 

26) C. Volpe, J. Phys. G34 (2007) Rl, |hep-ph/0605033| 

27) M. Lindroos, M. Mezzetto, London, UK: Imperial College Pr. (2010) 154 p. 

28) A. B. Balantekin, J. Gava and C. Volpe, Phys. Lett. B662 (2008) 396, 
larXiv:0710.31121 

29) J. Gava and C. Volpe, Phys. Rev. D78 (2008) 083007, larXiv:0807.34l8l 

30) J. Gava, C. Volpe, Nucl. Phys. B837 (2010) 50, larXiv:1002.098Tl 

31) D. Autiero, et al, JCAP 0711 (2007) 011. larXiv:0705.0116l 

32) D.Vaananen, C. Volpe, arXiv:1105.6225[ to appear in JCAP. 



33) C. Volpe, J. Phys. G30 (2004) LI, hep-ph/0303222 



